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Summary 

This  report  investigates  the  effects  of  dynamic  stall  and  three-dimensional  wake  on  the 
trim  results  of  control  inputs  and  periodic  responses  and  on  the  damping  levels  of  isolated 
hingeless  rotors.  The  investigation  covers  three  items:  1)  the  convergence  characteristics 
of  trim  and  damping  with  respect  to  the  number  of  harmonics  of  the  wake  model  in  the 
presence  of  dynamic  stall;  2)  a  parametric  study  of  trim  and  damping  over  a  broad  range 
of  system  parameters  such  as  thrust  level  and  advance  ratio,  and  3)  a  comprehensive  corre¬ 
lation  with  the  measured  lag-regressive-mode  damping.  The  correlation  includes  near-zero 
thrust  conditions  in  hover  and  forward  flight  to  high-thrust  and  highly  stalled  conditions 
in  forward  flight  with  advance  ratio  as  high  as  0.55  and  shaft  angle  as  high  as  20°.  The 
convergence  characteristics  and  the  parametric  study  are  based  on  a  three-bladed  isolated 
rotor  in  propulsive  trim;  the  cantilever  blades  have  flap  bending,  lag  bending  and  torsion 
degrees  of  freedom.  In  the  correlation  study,  the  root-flexure-blade  assembly  of  the  experi¬ 
mental  rotor  is  modeled  by  an  offset-hinged,  rigid  flap-lag  model  and  by  three  elastic  blade 
models,  which  differ  in  modeling  the  root  flexure.  The  experimental  rotor  is  soft  inplane  and 
operated  untrimmed.  The  ONERA  dynamic  stall  models  of  lift,  drag  and  pitching  moment 
and  a  finite-state  three-dimensional  wake  model  are  used.  The  control  inputs  and  periodic 
responses  as  well  as  the  Floquet  transition  matrix  about  that  response  are  obtained  by  pe¬ 
riodic  shooting  with  damped  Newton  iteration.  The  damping  levels  are  generated  from  a 
full  Floquet  analysis  that  includes  all  the  structural  and  aerodynamic  states  in  trim  analysis 
and  eigenanalysis.  For  lag  damping  predictions  in  the  presence  of  dynamic  stall,  at  least 
nine  wake  harmonics  are  required  even  for  the  low-frequency  lag-regressive  mode.  Moreover, 
dynamic  stall  and  wake  effects  are  appreciable  on  control  inputs,  periodic  responses  and 
damping  levels,  and  including  these  effects  appreciably  impi'oves  the  correlation. 
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Notation 


Unless  otherwise  stated,  the  symbols  below  are  dimensionless: 
a  linear  lift  curve  slope,  {rad~^) 

ad,  am  damping  factors  in  dynamic  stall  drag  and  pitching  moment  models 

b  airfoil  semi-chord,  (l/R) 

c,  c  airfoil  chord,  (m);  airfoil  chord,  (l/R) 

Cd,Cdo  airfoil  sectional  drag  coefficient  and  constant  profile  drag  coefficient 

Cl,  Cm  airfoil  lift  and  pitching  moment  coefficients 

Cmo  airfoil  pitching  moment  coefficient  at  zero  angle  of  attack 

Ct  thrust  coefficient 

Ct/cTs  thrust  level 

e  phase  shift  parameter  in  dynamic  stall  lift  model 

Ed,  Em  phase  shift  parameters  in  dynamic  stall  drag  and  pitching  moment  models 
/  equivalent  flat  plate  area  of  parasite  drag 

j  polynomial  number 

k  hjx 

Li  blade  sectional  lift  of  the  f-th  blade,  {Njm) 

Lx,Ly  local  aerodynamic  forces  (Fig.  1),  {Ij phLi^R?) 

Lq  apparent  mass  lift  normal  to  the  chord  line,  {Ij phVL^R?) 

Lc,  Lj  influence  coefficient  matrices 

m  harmonic  number 

M  aerodynamic  pitching  moment,  (l//)6n^i?‘*c);  also  number  of  harmonics 

in  the  wake  model 

Mq  noncirculatory  aerodynamic  pitching  moment  per  unit  length,  {\j phLRR^c) 

M  mass  matrix 

n  polynomial  number 

{ji)  normalized  Lengendre  polynomial  {n,m)  of  the  first  kind 
Q  number  of  blades 

r  harmonic  number 
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frequency  parameters  in  dynamic  stall  drag  and  pitching  moment  models 

radial  station  of  the  f-th  blade,  (l/R) 

rotor  radius,  {m) 

number  of  radial  shape  functions 

time  parameter 

resultant  air  velocity  at  a  blade  station  (Fig.  1),  {1/ClR) 
inflow  parameter 
total  inflow  parameter 

first  and  second  lag  mode  generalized  coordinates 
diag  {VuV,V,...,V} 

dynamic  stall  lift  frequency  and  damping  parameters 
first  and  second  flap  mode  generalized  coordinates 
radial  distance  measured  from  the  rotor  center,  {l/R) 
blade  airfoil  section  angle  of  attack,  {rad) 
wake  states 

shaft  tilt  angle,  {deg,  rad) 

Lock  number  (blade  inertia  parameter) 

circulation  per  unit  length,  UCi 

circulation-like  drag  per  unit  length,  UCd 

circulation-like  pitching  moment  per  unit  length,  UCm. 

pitch-rate  coefficient,  {rad~^ ) 

pitching-moment  parameter 

airfoil  rotation  rate  with  respect  to  airmass 

collective  pitch  angle,  {deg,  rad) 

longitudinal  cyclic  pitch  angle,  {rad) 

lateral  cyclic  pitch  angle,  {rad) 

time-delay  parameter 

downwash  on  the  i-th  blade 

thrust  induced  inflow 

total  inflow 
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advance  ratio 

first  and  second  torsion  mode  generalized  coordinates 
radial  shape  functions 

(n)  =  in  = 

damping  level  or  real  part  of  characteristic  exponent,  (1/sec) 
rotor  solidity 

azimuthal  location  of  the  f-th  blade 
air  density  (kg/m^) 
flap  natural  frequency,  rotating 
lag  natural  frequency,  rotating 
torsion  natural  frequency,  rotating 
rotor  angular  speed,  {rail sec) 
inflow  forcing  functions 
unstalled  and  stalled  components 
X  and  y  components 
time  derivative  of  (  ) 
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1  Introduction 


The  current  capability  to  predict  helicopter  lag  damping  needs  significant  improvements. 
Four  considerations  make  the  prediction  difficult  and  motivate  the  present  investigation  based 
on  Floquet  theory.  First,  drag,  induced  drag  and  Coriolis  forces  are  delicately  balanced  in  the 
inplane  direction.  Therefore,  from  hovering  to  high-speed  flights,  a  wake  representation  that 
goes  well  beyond  dynamic-inflow  and  lift-deficiency-type  approximations  is  required  (Ref.  1). 
Second,  with  increasing  demand  on  the  stability  margins  of  high-speed  and  highly  maneu¬ 
verable  helicopters,  high-thrust  conditions  in  forward  flight  become  increasingly  important. 
Thus,  the  complex  phenomenon  of  dynamic  stall  becomes  an  issue  (Refs.  2-4)  and  needs  to 
be  accounted  for  as  well.  Third,  lag-damping  prediction  is  sensitive  to  modeling  the  rotor 
system  and  its  flow  field.  It  is  also  sensitive  to  the  accuracy  of  computing  the  control  inputs 
and  the  corresponding  periodic  forced  responses  from  trim  analysis,  the  Floquet  transition 
matrix  about  that  trim  solution  and  the  eigenvalues  of  this  matrix.  Thus,  it  is  important 
that  a  consistent  sophistication  is  maintained  from  modeling  to  trim  analysis  to  eigenanal- 
ysis.  Often  this  requires  that  all  the  structural  and  aerodynamic  states  are  included  not 
only  in  trim  analysis  but  also  in  eigenanalysis  (Ref.  5).  This  is  a  demanding  and  expensive 
exercise;  with  lag-bending,  flap-bending  and  torsional  degrees  of  freedom,  even  simplified 
research  models  of  multiblade  rotors  with  wake  and  stall  dynamics  routinely  lead  to  nearly 
300  structural  and  aerodynamic  states  for  converged  results.  Fourth,  given  the  complexity 
of  the  rotor  system  and  its  flow  field,  it  becomes  necessary  that  the  predictions  are  checked 
out  against  a  broad  range  of  model  test  data.  Curi’ently  such  test  data  are  severely  limited. 
Although  not  typical  of  operational  rotors,  the  database  due  to  McNulty  is  an  exception 
(Ref.  6).  It  has  nearly  2000  data  points  and  is  comprehensive  to  include  near-zero  thrust 
conditions  in  hover  and  forward  flight  as  well  as  high-thrust  and  highly  stalled  conditions  in 
forward  flight  (0  <  ^  <  0.55,0°  <  0o  <  6°,0°  <  as  <  20°);  for  details  see  Refs.  2  and  7. 

Given  this  background,  we  review  the  aeroelastic  stability  studies  with  dynamic  stall 
and/or  three-dimensional  wake  modeling.  The  use  of  both  stall  and  three-dimensional  wake 
modeling  is  almost  routine  in  loads  and  vibration  studies.  By  comparison  their  use  is  far 
less  common  in  stability  studies.  Perhaps  the  first  study  to  consider  their  use  is  due  to 
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Torok  and  Chopra  (Ref.  4),  who  also  provide  a  review  of  aeroelastic  stability  modeling 
through  1991.  However,  in  Ref.  4,  stall  and  free  wake  are  included  in  trim  analysis,  but 
they  are  frozen  in  dynamic  perturbation.  Thus,  the  important  effects  of  wake  and  stall  on 
trim  analysis  are  included,  but  their  effects  on  eigenanalysis  are  not.  Moreover,  Ref.  4  uses 
only  a  lightly  loaded  section  of  the  database  of  Ref.  6  (shaft  tilt  <  10°  and  collective 
pitch  6q  <  3°).  Thus,  an  appreciable  section  of  the  database  under  high  shaft-angle  and 
collective-pitch  conditions  (a^  <  20°  and  Oq  <  6°)  is  not  included.  Developments  since  1991 
include  the  works  of  Barwey  et  al.  (Refs.  2  and  3),  who  use  dynamic  stall  modeling  and 
provide  correlation  with  virtually  the  complete  database  of  Ref.  6.  While  a  rigid  flap-lag 
model  is  used  in  Ref.  2,  two  elastic  flap-lag-torsion  models  that  differ  in  simulating  the  root 
flexure  are  used  in  Ref.  3.  Manjunath  et  al.  (Ref.  1)  use  a  finite-state  three-dimensional 
wake  model  and  rigid  flap-lag  modeling.  Reference  1  also  includes  limited  correlation  with 
test  data  of  Ref.  6  under  low-thrust  conditions.  The  work  in  Ref.  1  is  extended  in  Ref.  8 
to  include  elastic  blade  modeling.  References  1  and  8  also  include  a  review  of  aeroelastic 
stability  predictions  with  the  use  of  dynamic  wake  modeling. 

In  this  report,  we  investigate  the  effects  of  dynamic  stall  and  wake  on  the  trim  prediction 
of  control  settings  and  periodic  responses  as  well  as  on  the  damping  levels  of  isolated  hingeless 
rotors.  We  begin  with  a  convergence  study;  that  is,  we  investigate,  in  the  presence  of 
dynamic  stall,  the  degree  of  wake  modeling  necessary  for  converged  results  of  trim  and 
damping.  We  then  present  a  parametric  study  to  demonstrate  the  effects  of  dynamic  stall 
and  three-dimensional  wake  on  trim  and  stability.  We  emphasize  that  all  the  structural 
and  aerodynamic  states  are  included  in  trim  analysis  and  eigenanalysis.  This  is  followed 
by  a  comprehensive  correlation  with  the  complete  database  of  Ref.  6.  To  help  clarify  the 
aerodynamic  aspects  of  the  stability  problem,  predictions  based  on  the  linear,  quasisteady 
stall  and  dynamic  stall  theories  are  included  in  the  parametric  and  correlation  studies.  In  the 
correlation,  the  blade  dynamics  is  described  by  a  rigid  blade  model  as  well  as  by  three  elastic 
blade  models,  which  differ  in  modeling  the  root-flexure  of  the  experimental  rotor  (details  to 
follow).  The  intent  here  is  not  so  much  to  compare  one  structural  model  with  the  other; 
rather  it  is  to  help  clarify  further  the  aerodynamic  aspects  without  the  correlation  being 
masked  by  the  sensitivity  of  the  predictions  to  structural  modeling  approximations. 
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We  use  the  ONERA  dynamic  stall  models  of  lift,  drag  and  pitching  moment  (Refs.  9 
and  10).  The  unsteady  wake  is  described  by  a  finite-state  three-dimensional  wake  model 
(Ref.  11).  These  stall  and  wake  models  are  well  suited  to  finite-state  representation  and 
perturbations  about  a  trim  solution  and  thereby  to  Floquet  eigenanalysis.  This  is  noteworthy 
in  comparison  with  a  range  of  other  models  {t.g.  free  wake  model),  which,  though  with  proven 
usage  in  vibration  and  performance  analyses,  do  not  lend  themselves  well  to  eigenanalysis. 
Thus  in  summary,  the  trim  and  damping  predictions  are  based  on  a  broad  range  of  structural 
and  aerodynamic  representations,  and  the  correlation  covers  a  comprehensive  database  (Ref. 
6).  Such  a  work  from  wake  modeling  detail  to  parametric  study  to  correlation  should  provide 
a  useful  reference  and  demonstrate  the  strengths  and  weaknesses  of  predicting  lag  mode 
damping  in  the  presence  of  dynamic  stall. 

2  Modeling 

A  cantilever  blade  model  is  used  in  the  convergence  and  parametric  studies  of  trim 
and  damping.  The  correlation  study  requires  an  adequate  model  of  the  root-flexure-blade 
assembly  of  the  three-blade  experimental  rotor,  which  was  designed  to  represent  a  simple 
model  of  a  hingeless  rotor  with  spring  restrained  flap-lag  hinges.  While  the  mass  and  stiffness 
distributions  are  essentially  uniform  for  the  blade  portion,  they  are  highly  nonuniform  for  the 
root  flexure  or  root  beam.  According  to  Ref.  3,  the  sensitivity  of  lag-damping  predictions 
to  the  structural  refinements  in  modeling  the  root-flexure-blade  assembly  increases  with 
increasing  blade  pitch  and  advance  ratio.  Therefore,  to  bracket  the  extent  of  this  sensitivity 
to  modeling  details,  we  use  four  models  of  the  root-flexure- blade  assembly.  They  range  from 
a  rigid  flap-lag  blade  model  to  three  elastic  flap-lag-torsion  models,  which  have  identical 
blade  representation  but  differ  in  root-flexure  representations.  All  four  models  have  the 
capacity  to  simulate  full  and  zero  structural  flap-lag  couplings  of  the  experimental  model. 
The  lag-damping  predictions  are  based  on  four  aerodynamic  theories,  which  range  from 
linear  (quasisteady)  theory  to  a  relatively  comprehensive  dynamic  stall  and  wake  theory.  To 
simplify  the  subsequent  presentation  of  results  and  clarify  terminology,  we  briefly  describe 
these  four  models  of  root-flexure-blade  assembly  and  four  aerodynamic  theories. 
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2.1 


Root-Flexure-Blade  Assembly 

2.1.1  Rigid  Flap-Lag  Model 

It  is  an  offset-hinged  rigid  blade  with  flap  and  lag  (inplane)  degrees  of  freedom;  the  flap 
and  lag  hinges  are  spring  restrained  and  coincident  (Ref.  2). 

2.1.2  Spring  Model 

The  root-flexure  is  simulated  by  a  set  of  three  linear  springs  located  at  an  effective 
hinge  offset;  the  spring  stiffnesses  are  based  on  measured  values.  The  airfoil  or  the  blade 
portion  has  uniform  mass  and  stiffness  distributions.  The  bending-torsion  couplings  of  the 
root  spring  system  are  fully  accounted  for  (Ref.  3). 

2.1.3  Modified  Model 

The  root  flexure  is  simulated  by  a  short  beam,  over  which  the  mass  and  stiffness  prop¬ 
erties  are  kept  uniform.  Their  magnitudes  are  chosen  so  as  to  reasonably  match  the  funda¬ 
mental,  nonrotating  bending  and  torsion  frequencies  of  the  experimental  model,  and  thus  to 
those  of  the  spring  model.  With  judicious  choice  of  hinge  offset  and  the  length  of  root  beam 
it  is  possible  to  match  the  nonrotating  and  rotating  frequencies  of  the  first  two  flap-bending, 
lag-bending  and  torsional  modes  of  the  spring  model  (Ref.  3). 

2.1.4  Modal  Model 

The  actual  mass  and  stiffness  distributions  of  the  root-flexure-blade  assembly  of  the  test 
model  are  used  in  calculating  the  nonrotating  mode  shapes  numerically. 

2.2  Aerodynamics 

The  aerodynamic  representation  called  the  dynamic  stall  and  wake  theory  is  fairly  com¬ 
prehensive.  It  includes  the  effects  of  dynamic  stall  lift,  drag  and  pitching  moment  from  a  thin 
airfoil  theory  (Refs.  9  and  10)  and  the  downwash  effects  from  a  finite-state  three-dimensional 
wake  theory  (Ref.  11).  Moreover,  the  airfoil  theory  includes  the  effects  of  reversed  flow  and 
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large  angles  of  attack,  and  the  wake  theory  accounts  for  the  finite  number  of  blades.  We 
begin  with  a  brief  account  of  dynamic  stall  theory.  This  is  followed  by  a  mention  of  how 
quasisteady  stall  and  linear  theories  are  recovered  as  special  cases.  We  conclude  this  section 
with  a  discussion  of  dynamic  stall  and  wake  theory. 

2.2.1  Dynamic  Stall  Theory 

It  is  based  on  the  ONERA  models  of  unified  lift,  drag  and  pitching  moment.  Basically 
we  introduce  lift  circulation  F,  circulation-like  drag  Fj  and  circulation- like  pitching  moment 
Fm  as  follows: 

r  =  ri-fr2,  =  Fm  =  Fm,  +  ^m-2 

where  subscripts  1  and  2  indicate  the  linear  and  stalled  or  nonlinear  components,  respectively. 
These  six  components,  two  each  in  lift,  drag  and  pitching  moment  are  governed  by  Eqs.  (2)  - 

(4). 

Dynamic  Stall  Lift: 

fcfi -H  AFi  =  Aa(C/y -1- 6£)cosa  +  (2a) 

-t-  2dwkt2  4-  (l-f  r2  =  (l  +  d^)  [U  +  ek  (tl^  cos  a  +  Uy  sin  a)  AC^ 


4-  ek  [Uy  sin  a  -  tlx  cos  a'j  ~~q^] 

Dynamic  Stall  Drag: 

'  Fd,  =  UCi,  (3a) 

-h  aikido  rlVd,  =  -  [rlU^Cd  -H  E^kUy]  (3b) 

Dynamic  Stall  Pitching  Moment: 

(4^) 

k^f^,  +  amkVm,  4-  rlVm,  =  -  [r^RAC,„  -f  Err^kUy]  (4b) 


It  is  seen  that  the  linear  components  in  Eqs.  (2a),  (3a)  and  (4a)  follow  the  classical  thin 
airfoil  theory.  By  comparison,  the  stalled  components  have  an  involved  algebraic  structure 
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and  merit  additional  comments.  For  example,  in  Eqs.  (2b),  (3b)  and  (4b),  ACz,ACd  and 
ACm  act  like  driving  forces  and  represent  the  differences  between  the  linear  and  quasisteady 
values  of  the  airfoil  section  characteristics.  For  example,  ACz  at  an  instantaneous  angle  of 
attack  a  is  the  difference  between  the  extrapolated  linear  lift  coefficient  and  the  quasisteady 
stall  lift  coefficient  of  the  airfoil;  for  details,  see  Ref.  12.  Similarly,  A,  6,  d,  e  and  w  in  the  lift 
equation;  ad,rd  and  Ed  in  the  drag  equation  and  and  Em  in  the  pitching-moment 

equation  are  determined  on  the  basis  of  wind-tunnel  experiments.  Another  important  pa¬ 
rameter  in  the  linear  part  of  the  lift  equation  is  i,  which  represents  the  airfoil  rotation  rate 
relative  to  airmass  and  includes  complete  geometric  rotations  of  the  airfoil. 

It  is  expedient  to  represent  the  above  airfoil  lift  F,  drag  and  pitching  moment  F^  in 
the  local  airfoil  coordinates  in  terms  of  Lx-,Ly  and  M  as  shown  in  Fig.  1.  This  is  done  in 
Eqs.  (5a)-(5c). 


Ty  =  f/j;  [Fl  -(-  r2]  -f  F’y  [r^j  -f  F(i2  ]  +  Fq 

(5a) 

L^  =  -Uy[Ti  +  T2]  +  UATd,+rd,] 

(5b) 

M  =  2b[U(Tm,+  Tm2)]  +  Mo 

(5c) 

where  Lq  and  Mq  are  apparent  mass  lift  normal  to  the  chord  and  noncirculatory  pitching 
moment  at  the  three-quarter  chord  point. 

2.2.2  Quasisteady  Stall  Theory 

It  includes  the  airfoil-section  quasisteady  stall  characteristics.  From  Eqs.  (2)-(4),  by 
suppressing  dynamic  stall  characteristics  we  get 

Fi  =  a  (f7y  4- 6^)  cos  a;  V2  =  -UACz  (6a) 


=  UCd,-,  Td,  =  -UACd 


(6b) 


2.2,3  Linear  Theory 


Similarly,  by  suppressing  the  quasisteady  stall  characteristics  the  equations  of  lift,  drag 
and  pitching  moment,  including  effects  of  reversed  flow  and  large  angles  of  attack,  are: 


T 1  a  {Uy  +  bi)  cos  a;  r2  =  0 

(7a) 

Fdj  =  UCdo',  =  0 

(7b) 

Tmi  —  UCrriQ  -|-  Fttij  —  0 

(7c) 

2.2.4  Dynamic  Stall  and  Wake  Theory 

It  is  dynamic  stall  theory  with  downwash  dynamics  modeled  by  a  finite-state  three- 
dimensional  wake  theory  due  to  Peters,  Boyd  and  He  (Ref.  11).  At  a  blade  station  with 
radial  coordinate  fj-  and  spatial  azimuth  xpi,  the  instantaneous  wake  or  downwash  A  (f;,  rl^i,t) 
is  given  by  a  complete  set  of  radial  shape  functions  (f,)  and  spatial  harmonics  cos{r'il^i) 
and  sin(r'0j): 

oo  oo 

A  (r,-,  t)  =  Y^  ^  <!>]  (r,)  [a5(t)  cos{ripi)  -f  sin(ri/)i)]  (8) 

r=0  j=r+l,r+3 

The  cosine  component  and  the  sine  component  are  the  dynamic  states  of  the 

downwash  and  are  governed  by 


M{dj}-fVL-i{a;:}  =  0.5{rr}  (9a) 

M  {/jj}  +  VL;^  {/Ij}  =  0.5  {r^}  (9b) 


where  V  is  the  diagonal  matri.x  with  Vii  =  =  y^(/i^  +  A^ )  and  all  other  elements  are 

given  by  y  =  [/^^  +  (Aj  -f-  A^)  At]  j +  Af).  Closed-form  expressions  are  available  for  the 
diagonal  mass  matrix  M  and  influence  coefficient  matrices  Lc  and  L5.  Similarly,  r™'"  and 
r™*  are  cosine  and  sine  components  of  the  pressure  coefficient,  which,  for  a  rotor  with  Q 
blades,  are  given  by 


-Oc 


(10a) 


j^mc 


pQ.'^R^ 


ch'i  cos 


(10b) 


i 


(10c) 


1  Q  /•] 

C'  =  iE/ 

TT^JO 


pQ.'^R? 


di'i  sin  {mtpi) 


A  final  comment  concerns  the  computation  of  equilibrium-state  infiow.  In  linear,  quasis¬ 
teady  stall,  and  dynamic  stall  theories  it  is  computed  from  the  momentum  theory,  while  in 
the  dynamic  stall  and  wake  theory  it  is  computed  from  the  coupled  blade-wake-stall  equa¬ 
tions. 


3  Analysis 

For  the  rigid  flap-lag  model  of  the  root-flexure-blade  assembly,  we  follow  Ref.  2  for 
the  equations  of  motion  including  the  hinge-offset  effects.  For  the  elastic  blade  models,  the 
equations  of  flap  bending,  lag  bending  and  torsion  are  based  on  Hamilton’s  principle  with  a 
second-order  ordering  scheme.  The  spatial  dependence  is  treated  by  a  Galerkin  scheme  us¬ 
ing  the  uncoupled  nonrotating  mode  shapes.  Development  of  closed-form  expressions  of  the 
mode  shapes  is  a  routine  exercise  for  the  cantilever-beam  and  spring  models.  However,  for 
the  modified  model,  we  use  the  closed-form  expressions  developed  in  Ref.  3  using  computer 
algebra;  the  mode  shapes  refer  to  a  stepped  beam  with  one  end  fixed  and  the  other  free.  At 
the  junction  of  the  blade  and  the  root  flexure,  there  is  discontinuity  in  the  distributions  of 
mass  and  flap-lag-torsion  stiffnesses.  The  closed-form  expressions  account  for  this  discon¬ 
tinuity  and  provide  continuous  displacements,  slopes,  moments  and  shear  forces.  For  the 
modal  model,  the  nonrotating  modes  are  calculated  by  Myklestad-type  and  finite  element 
schemes  with  mass  and  stiffness  distributions  identical  to  those  of  the  e.xperimental  model. 
The  modes  from  these  two  schemes  agree.  Moreover,  we  have  also  done  extensive  spot  checks 
using  the  rotating  mode  shapes;  trim  and  damping  predictions  from  these  two  types  of  mode 
shapes  agree.  Throughout,  the  Galerkin- type  integrals  are  calculated  numerically. 

As  in  the  experimental  rotor  (Refs.  6  and  7),  the  analytical  model  in  the  correlation  has 
three  blades  of  the  NACA  23012  airfoil  .section  and,  as  detailed  earlier,  the  root-flexure-blade 
assembly  modeling  varies  from  a  simple  offset-hinged  rigid  blade  model  to  three  different 
root-flexure  models  with  elastic  blades.  In  the  convergence  study  on  wake  modeling  and  in 
the  parametric  study  we  consider  an  isolated  rotor  with  three  cantilever  blades  of  the  NACA 
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23012  airfoil  section.  To  include  dynamic  stall  lift,  drag  and  pitching  moment,  the  blade  is 
discretized  into  five  equal-length  elements;  each  element  has  7  dynamic  stall  states  (3  in  lift 
and  2  each  in  drag  and  pitching  moment).  The  airfoil-section  quasisteady  stall  and  dynamic 
stall  characteristics  of  lift,  drag  and  pitching  moment  for  a  complete  sweep  of  angle  of  attack 
(—180°  <  cx  <  180°)  are  as  in  Ref.  3.  In  the  convergence  and  parametric  studies,  the  trim 
analysis  refers  to  computing  the  control  inputs  —  collective  pitch  angle  Oq,  longitudinal  cyclic 
pitch  angle  9ci  lateral  cyclic  pitch  angle  0^  and  shaft  tilt  angle  ag  —  and  the  corresponding 
periodic  responses.  The  control  inputs  and  the  initial  conditions  for  periodic  response  are 
computed  simultaneously  by  periodic  shooting.  The  resulting  nonlinear  equations  are  solved 
by  a  damped  Newton  iteration  with  a  damping  parameter  of  0.7.  Occasional  divergence 
problems  are  encountered  at  high-thrust  and  high-advance-ratio  conditions.  Under  such 
conditions,  the  results  are  obtained  by  either  adjusting  the  Newton  damping  parameter  or 
using  a  continuation  approach  with  no  damping;  also  see  Ref.  13.  The  Floquet  transition 
matrix  for  perturbations  about  periodic  response  comes  out  as  a  byproduct  of  periodic 
shooting,  and  its  eigenanalysis  gives  the  damping  levels.  All  the  structural  and  aerodynamic 
states  are  included  in  trim  analysis  as  well  as  in  eigenanalysis.  In  the  correlation  study, 
the  collective  pitch  angle  9q  and  shaft  tilt  angle  are  known  control  inputs  and  the  trim 
analysis  refers  to  computing  the  periodic  responses. 

4  Results 

The  results  are  presented  in  three  parts:  i)  wake-modeling  detail  necessary  for  converged 
results  of  trim  and  damping;  ii)  parametric  study  of  these  results  to  variations  in  system 
parameters  and  in)  correlation  of  damping  predictions  with  the  database  of  Ref.  6.  In  the 
first  two  parts,  we  consider  the  rotor  in  propulsive  or  flight  trim  and  include  all  multiblade 
lag  modes.  In  the  third  part,  the  correlation  covers  nearly  the  entire  database  of  some  2000 
points  on  lag  regressive-mode  damping.  Unless  otherwise  specified,  the  following  baseline 
values  are  used  in  the  first  two  parts:  Q  —  'i]  cr  —  0.05;  a  =  5.73;  7  =  5.29;  =  1.15;  = 

0.7;  u;0  =  5.0;  Cd^  =  0.0079;  /  =  0.01.  In  the  third  part,  the  test  focuses  on  the  forward-flight 
aspects  of  the  stability  problem.  The  rotor  is  operated  untrimmed  with  an  unrestricted  tilt 
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of  the  tip-path  plane.  To  cover  the  gist  of  the  database,  it  becomes  necessary  to  include 
samples  of  all  typical  cases  tested.  Accordingly,  we  include  the  following  cases:  collective 
pitch  angle  (0°  <  <  6°),  shaft  angle  (8°  <  tij  <  20°)  and  advance  ratio  {0  <  n  <  0.55). 

The  rotor  speed  Q,  =  1000  rpm,  corresponding  to  a  rotating  lead-lag  frequency  =  0.61 
and  flap  frequency  up  =  1.1.  The  flap-lag  structural  coupling  is  zero;  this  refers  to  the 
database  for  which  the  collective  pitch  angle  is  set  manually  by  changes  to  the  angle  of  the 
blade  relative  to  the  root  flexure.  The  thrust  level  is  obtained  as  part  of  the  trim  analysis. 
Low-thrust  conditions  occur  when  both  the  collective  pitch  angle  and  advance  ratio  are  low 
or  when  a  high  collective  pitch  angle  offsets  a  high  shaft  angle.  Similarly,  high  (negative) 
thrust  conditions  occur  at  a  low  collective  pitch  angle  when  both  the  advance  ratio  and  shaft 
angle  are  high.  Thus,  the  correlation  covers  near-zero-  to  high-thrust  conditions  for  various 
combinations  of  shaft  angle,  collective  pitch  angle  and  advance  ratio.  For  additional  details 
see  Refs.  2,  6  and  7. 

4.1  Wake  Modeling 

We  now  come  to  wake  modeling  for  converged  results  of  control  inputs,  periodic  responses 
and  damping  levels.  This  involves  a  judicious  combination  of  the  number  of  radial  shape 
functions  with  each  harmonic  or  S  and  the  total  number  of  harmonics  or  M;  see  Eq.  (8). 
It  is  less  sensitive  to  S  per  se,  M  being  the  more  dominant  parameter.  Results  (not  shown) 
demonstrate  that  three  radial  shape  functions  with  each  harmonic  ai’e  adequate,  as  was  the 
case  in  Refs.  1  and  8,  which  are  based  on  linear  aerodynamics;  that  is,  linear  quasisteady 
theory  with  wake  dynamics.  Accordingly,  with  S  =  3,  we  determine  M.  We  compute  the 
control  inputs,  periodic  responses  and  corresponding  damping  levels  for  a  preset  value  of  M 
and  then  we  increase  M  in  increments  of  one  till  convergence  to  preset  tolerance  is  obtained. 
For  Cr/cTs  =  0.075  and  fx  =  0.3,  Figs.  2,  and  3  show,  respectively,  the  control  inputs,  and 
the  damping  levels.  In  particular.  Fig.  2  shows  all  four  control  inputs  —  three  blade  pitch 
settings  and  shaft  tilt.  It  is  seen  that  M  =  9  essentially  gives  converged  trim  results  in 
that  maximum  error  with  respect  to  the  results  for  M  =  12  is  less  than  4%.  Similarly,  Fig. 
3  shows  the  damping  levels  of  all  three  multiblade  lag  modes.  We  observe  an  oscillatory 
convergence  trend,  particularly  with  respect  to  the  cyclic  regressive  and  progressive  modes, 
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and  M  =  9  gives  the  same  degree  of  convergence  as  that  obtained  for  the  control  inputs  in 
Fig.  2.  In  Figs.  2  and  3,  the  case  with  M  =  I  merits  special  mention;  it  is  a  dynamic- inflow 
approximation  and  the  results  clearly  show  the  inadequacy  of  this  approximation.  Similarly, 
convergence  of  control  inputs,  periodic  responses  and  the  corresponding  damping  levels  are 
investigated  for  different  values  of  Cr/cr^  and  /i  covering  the  flight  regimes  of  conventional 
helicopters.  Although  convergence  degrades  with  increasing  and  /i,  the  wake  model 

with  5  =  3  and  M  =  9  (57  wake  states)  is  adequate. 

4.2  Parametric  Study 

Figures  4-10  show  parametric-study  results  with  an  essentially  converged  wake  model 
(57  wake  states  with  M  =  9  and  S  =  3).  The  results  are  based  on  four  aerodynamic 
theories,  which  are  identified  in  the  figures.  Figure  4  shows  the  periodic  responses  of  all 
six  modal  components  —  two  each  in  flap  bending,  lag  bending  and  torsion;  advance  ratio 
fi  =  0.3  and  thrust  level  CtIcTs  =  0.1.  Appreciable  quantitative  differences  are  seen  among 
these  theories.  For  example,  the  amplitudes  of  the  first-lag  (first-flap)  modes  are  0.00324 
(0.00505),  0.00548  (0.006),  and  0.00646  (0.00775)  from  the  linear,  quasisteady  stall  and 
dynamic  stall  theories,  respectively.  Compared  to  0.009  (0.0076)  from  the  dynamic  stall  and 
wake  theory,  the  respective  percentage  errors  are  64  (34),  39  (21)  and  28  (2).  The  linear 
theory  vis-a-vis  the  quasisteady  stall  theory  shows  the  effects  of  nonlinear  airfoil-section 
characteristics.  Similarly,  we  observe  appreciable  effects  of  dynamic  stall  and  dynamic  wake 
by  comparing,  respectively,  dynamic  stall  theory  with  the  quasisteady  stall  theory  and  with 
the  dynamic  stall  and  wake  theory.  Figure  4  also  shows  that  dynamic  stall  and  dynamic 
wake  bring  in  qualitative  changes  in  the  periodic  responses  of  high-frequency  (>  2/rev) 
modes.  For  example,  dynamic  stall  theory  and  dynamic  stall  and  wake  theory  predict  4/rev 
oscillations  in  the  second  flap  mode  and  5/rev  oscillations  in  the  first  torsion  mode.  Thus, 
Fig.  4  demonstrates  the  necessity  of  including  the  effects  of  dynamic  stall  and  wake  in  the 
prediction  of  periodic  responses  in  trim  analysis.  These  aerodynamic  aspects  are  further 
investigated  in  Figs.  5-7  with  respect  to  the  four  control  inputs.  The  thrust  level  Ct/ct^ 
is  0.075  in  Fig.  5  and  0.1  in  Fig.  6  and  the  advance  ratio  n  varies  from  hover  to  0.4  in 
both  the  figures.  Figure  7  is  a  cross-plot  of  these  two  figures  for  varying  thrust  level  at 
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pL  =  0.4.  Overall,  these  three  figures  show  that  for  jj,  <  0.4  and  at  low-thrust  conditions 
{CtIcTs  <  0.05)  all  four  theories  are  fairly  close;  this  is  expected.  However,  with  increasing 
thrust  level  and  advance  ratio  the  dynamic  stall  and  dynamic  wake  effects  increase,  and 
so  do  the  differences  among  these  theories.  For  example,  at  fi  =  0.4  and  Cx/cr^  =  0.1 
the  longitudinal  cyclic  pitch  6c  is  equal  to  3°  (or  0.05  rad)  and  9°  according  to  linear  and 
quasisteady  stall  theories  whereas  it  rapidly  increases  to  23°  and  27°  according  to  dynamic 
stall  and  dynamic  stall  and  wake  theories.  This  is  because  a  major  portion  of  the  rotor  disk 
is  stalled  under  such  high-thrust  and  high-advance  ratio  conditions;  this  demands  high  pitch 
settings  to  maintain  the  required  thrust  and  moment  conditions.  Moreover,  compared  to  the 
sensitivity  of  three  pitch  settings  to  dynamic  stall  and  dynamic  wake,  the  shaft  tilt  is  less 
sensitive.  Overall,  Figs.  5-7  show  that  the  effects  of  dynamic  stall  and  wake  rapidly  increase 
with  increasing  thrust  level  and  advance  ratio. 

Figures  8-10  show  the  damping  levels  of  the  three  multiblade  lag  modes;  the  three- 
bladed  rotor  is  in  propulsive  trim  and  the  predictions  are  based  on  the  four  aerodynamic 
theories  referred  to  earlier.  The  thrust  level  Cx/cTj  =  0.075  in  Fig.  8  and  it  increases 
to  0.1  in  Fig.  9.  Figure  10  is  a  sort  of  cross-plot  of  Figs.  8  and  9  at  a  relatively  high 
advance  ratio  of  0.4,  and  it  includes  additional  damping  results  for  a  thrust  level  of  0.05.  As 
seen  from  Fig.  8,  the  quasisteady  stall  theory  deviates  appreciably  from  the  linear  theory; 
by  contrast  it  is  closer  to  the  other  two  theories:  dynamic  stall,  and  dynamic  stall  and 
wake  theories.  This  deviation  is  due  to  nonlinear  airfoil  characteristics  of  substall  drag 
and  quasisteady  stall  lift  and  drag.  The  effects  of  quasisteady  stall  pitching  moment  are 
negligible.  At  low  advance  ratios,  the  substall  drag  dominates,  and  with  increasing  advance 
ratio  the  effects  of  quasisteady  stall  lift  and  drag  become  increasingly  important  in  the 
quasisteady  stall  theory.  Up  to  an  advance  ratio  of  0.25  or  so,  the  quasisteady  stall  theory 
agrees  with  the  dynamic  stall  theory  and  to  some  extent  with  the  dynamic  wake  and  stall 
theory  as  well.  From  hover  to  an  advance  ratio  of  0.1,  dynamic  stall  and  wake  theory 
predicts  less  damping  for  the  regressive  and  progressive  modes  owing  to  downwash  effects. 
For  advance  ratio  fi  >  0.25,  the  quasisteady  stall  theory  consistently  underpredicts  damping, 
and  this  underprediction  increases  with  increasing  advance  ratio  owing  to  increasing  dynamic 
stall  effects;  the  minor  differences  between  the  dynamic  stall  and  dynamic  stall  and  wake 
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theories  are  due  to  downwash  dynamics.  Figure  9  is  basically  a  magnified  version  of  Fig.  8, 
particularly  regarding  the  effects  of  dynamic  stall  and  wake  with  increasing  thrust  level.  The 
regressive  mode  damping,  in  particular,  shows  that  dynamic  wake  effects  are  appreciable  for 
practically  all  advance  ratios  {0  <  /.i  <  0.4).  A  comparison  between  the  quasisteady  stall 
and  dynamic  stall  theories  shows  that  for  /i  >  0.2  dynamic  stall  effects  rapidly  increase  with 
increasing  advance  ratio.  Figure  10  is  more  or  less  a  summary  of  Figs.  8  and  9,  although 
for  a  fixed  value  of  fi.  Even  at  low-thrust  levels  <  0.05)  the  inadequacy  of  the  linear 

theory  in  predicting  lag  damping  is  clearly  seen  and  as  expected,  the  other  three  theories 
agree  well.  By  comparing  dynamic  stall  theory,  respectively,  with  the  quasisteady  stall  theory 
and  with  the  dynamic  stall  and  wake  theory,  Fig.  10  shows  that  the  effects  of  dynamic  stall 
and  wake  increase  with  increasing  thrust  level.  Overall,  Figs.  8-10  demonstrate  the  necessity 
of  including  dynamic  stall  and  wake  in  predicting  damping  levels. 

4.3  Correlation 

We  reiterate  that  the  damping  levels  are  generated  from  four  structural  models  of  the 
root-flexure-blade  assembly  and  that  the  predictions  from  each  structural  model  are  based 
on  the  four  aerodynamic  theories  referred  to  earlier.  With  such  a  wide  range  of  structural 
and  aerodynamic  representations,  we  bracket  the  extent  of  the  sensitivity  of  the  predictions 
to  structural  approximations  and  thereby  isolate  the  aerodynamic  aspects  of  the  delicate  lag¬ 
damping  predictions  without  allowing  the  correlation  to  be  overly  masked  by  this  sensitivity. 
The  computer  time  for  such  comprehensive  predictions  covering  nearly  2000  test  points  of 
the  database  on  lag-regressive  mode  damping  is  nearly  prohibitive.  To  somewhat  minimize 
this  constraint,  the  correlation  uses  a  wake  model  with  5"  =  3  and  M  =  1  (not  9  as  done 
in  Figs.  2-10)  and  the  maximum  error  compared  to  the  results  with  S  =  Z  and  M  =  9 
is  less  than  12%.  In  Figs.  11-13  we  present  the  correlation  based  on  the  modified. model, 
and  then  in  Figs.  14-16  we  compare  the  predictions  from  the  four  structural  models  of  the 
root-flexure-blade  assembly  using  the  dynamic  stall  and  wake  theory. 

Figures  11a— lid  show  correlation  for  Oq  =  0°  at  high  shaft  angles;  as  =  12°,  14°,  16° 
and  20°,  respectively.  As  seen  from  the  data,  the  damping  remains  nearly  constant  from 
hover  to  an  advance  ratio  of  /z  ft;  0.2  and  then  increases  sharply  beyond  yu  =  0.3.  From 
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hover  to  an  advance  ratio  of  0.1,  owing  to  very  low-thrust  conditions,  the  predictions  from 
all  four  aerodynamic  theories  are  virtually  identical.  For  ^  >  0.3,  some  of  the  theories  differ 
qualitatively.  This  result  is  expected;  for  the  present  combination  of  low  collective  pitch  and 
high  shaft  tilt  angles,  as  the  advance  ratio  increases,  (negative)  thrust  level  increases  and 
thereby  dynamic  stall  dominates.  Basically,  the  linear  theory  predicts  that  damping  increases 
slowly,  and  thereafter  it  significantly  deviates  from  the  data,  leveling  off  and  decreasing 
sharply.  It  is  this  sharp  decrease  that  is  not  supported  by  the  data.  While  this  decrease 
occurs  within  the  data-range  for  shaft  angles  12°  and  14°,  it  occurs  outside  the  data-range 
for  shaft  angles  16°  and  20°.  The  quasisteady  stall  theory  shows  sharply  decreasing  damping 
for  0.35  <  fi  <  0.45  at  =  12°  and  14°,  and  for  0.3  <  /i  <  0.4  at  =  16°  and  20°,  and 
thereafter  it  shows  rapidly  increasing  damping.  The  increase  is  so  delayed  that  it  makes  the 
quasisteady  stall  theory  unacceptable.  In  contrast,  we  have  the  dynamic  stall,  and  dynamic 
stall  and  wake  theories.  For  example,  inclusion  of  dynamic  stall  dramatically  improves  the 
correlation.  The  dynamic  stall  theory  follows  the  quasisteady  stall  theory  till  n  w  0.3,  and 
for  >  0.3  the  dynamic  stall  effects  dominate.  In  other  words,  for  //  >  0.3,  the  dynamic  stall 
theory  deviates  remarkably  from  the  quasisteady  stall  theory  and  improves  the  correlation. 
The  inclusion  of  three-dimensional  wake  effects  in  addition  to  stall  effects  brings  in  further 
improvement.  Overall,  the  dynamic  stall,  and  dynamic  stall  and  wake  theories  predict  the 
trend  of  the  data;  that  is,  the  damping  basically  increases  with  increasing  advance  ratio, 
although  it  is  accompanied  by  localized  humps  and  troughs. 

Figure  12  shows  the  correlation  for  Oq  =  3°.  In  Fig.  12a,  =  8°  and  in  Fig.  12b, 

as  =  16°;  the  data  are  available  for  0  <  /f  <  0.25  and  for  0  <  <  0.35,  respectively. 

In  Fig.  12a,  the  thrust  level  is  relatively  low  throughout  (Cj/a^  <  0.04).  It  is  about 
0.03  in  hover,  and  increases  to  0.04  at  /«  =  0.1  and  thereafter  continuously  decreases  with 
increasing  advance  ratio,  down  to  nearly  zero  at  /4  w  0.55.  This  is  well  reflected  by  the 
good  correlation  of  the  linear  theory  and  closeness  of  the  quasisteady  stall  and  dynamic  stall 
theories.  Though  the  quasisteady  stall  and  dynamic  stall  theories  overpredict  damping,  they 
predict  the  trend  of  the  data  of  nearly  constant  damping  for  the  entire  advance  ratio  range 
(0  <  ^  <  0.25).  The  dynamic  stall  and  wake  theory  falls  in  between  the  linear  theory,  and 
the  quasisteady  stall  and  dynamic  stall  theories.  It  improves  the  correlation  by  significantly 
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reducing  the  overprediction  of  the  quasisteady  stall  and  dynamic  stall  theories  and  thus 
provides  satisfactory  correlation  throughout. 

The  data  in  Fig.  12b  show  clearly  the  trend  of  decreasing  damping  with  increasing 
advance  ratio.  The  thrust  level,  which  is  about  0.03  in  hover,  changes  to  —0.05  at  =  0.35. 
The  linear  theory  predicts  that  the  damping  consistently  increases  with  increasing  advance 
ratio,  a  trend  opposite  to  that  of  the  data.  However,  for  >  0.5  it  predicts  the  trend  of 
decreasing  damping  with  increasing  advance  ratio  ;  but  this  decrease  occurs  far  beyond  the 
data  range  to  be  adequate.  The  quasisteady  stall  and  dynamic  stall  theories  come  close 
to  predicting  the  trend  of  the  data  for  0  <  fi  <  0.25,  although  the  damping  levels  are 
overpredicted  owing  to  nonlinear  substall  drag  effects.  The  dynamic  stall  and  wake  theory 
shows  improvement  over  the  other  three  theories;  still  it  falls  short  of  predicting  clearly 
the  trend  of  the  data  within  the  data  range.  Given  the  very  low  thrust  condition,  neither 
dynamic  stall  nor  dynamic  wake  is  a  major  factor,  and  the  required  improvements  merit 
further  study. 

Now,  we  present  the  correlation  for  Oq  =  6°  and  =  20°  in  Fig.  13.  The  linear 
theory  shows  that  damping  increases  with  increasing  advance  ratio  from  hover  to  n  =  0.45 
followed  by  a  sharp  decrease  for  /i  >  0.5.  This  is  opposite  to  the  trend  of  the  data,  which, 
though  limited  to  very  low  advance  ratios  (0  <  /i  <  0.15),  show  that  damping  very  slowly 
decreases  with  increasing  advance  ratio.  The  quasisteady  stall  and  dynamic  stall  theories 
predict  the  trend  of  the  data  despite  appreciable  overpredictions.  The  dynamic  stall  and 
wake  theory  significantly  reduces  these  overpredictions  and  provides  good  correlation  over 
the  data  range.  The  nonlinear  substall  drag  accounts  for  much  of  the  differences  between  the 
linear  and  quasisteady  stall  theories;  this  is  a  consequence  of  low- Reynolds- number  effects 
of  the  test  conditions.  Dynamic  stall  is  not  an  issue  here,  and  unsteady  lift  effects  account 
for  the  differences  between  the  dynamic  stall  and  quasisteady  stall  theories.  Similarly,  the 
effects  of  downwash  account  for  the  differences  between  the  dynamic  stall,  and  dynamic  stall 
and  wake  theories. 
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4.3.1  Comparison  of  Structural  Models 

The  preceding  correlations  based  on  the  modified  model  (Figs.  11  —  13)  and  additional 
correlations  based  on  other  structural  models  (see  Ref.  14)  of  the  root-flexure-blade  as¬ 
sembly  show  that  the  dynamic  stall  and  wake  theory  is  the  best  of  the  four  aerodynamic 
theories.  Given  this  background,  we  present  in  Figs.  14—16  a  comparative  assessment  of 
four  structural  models  based  on  the  dynamic  stall  and  wake  theory.  Figures  14a— 14d  show 
the  effects  of  structural  modeling  for  Oq  =  0*^  and  =  12°,  14°,  16°  and  20°.  To  set  the 
stage,  we  mention  that  with  increasing  advance  ratio  and  shaft  angle,  the  negative  thrust 
level  increases,  as  do  the  effects  of  dynamic  stall.  This  brings  in  increased  system  nonlin¬ 
earity,  and  accordingly  damping  predictions  show  increasing  sensitivity  to  changes  in  blade 
modeling.  This  sensitivity  is  more  clearly  seen  in  Figs.  14c  and  14d,  for  a*  =  16°  and  20°, 
respectively.  For  shaft  angles  12°  and  14°,  the  predictions  from  three  elastic  models  are  close 
and  provide  adequate  correlation  throughout.  The  rigid  flap-lag  model  also  provides  satis¬ 
factory  correlation  up  to  an  advance  ratio  of  /i  0.3.  The  difference  between  the  predictions 
from  three  elastic  blade  models  and  those  from  the  rigid  blade  model  is  due  to  blade  torsion 
and  bending  effects  (Ref.  3).  Overall,  Figs.  14a-14d  show  that  all  the  elastic  blade  models 
adequately  predict  the  trend  of  the  data. 

Figure  15a  for  9q  =  3°  and  =  8°  shows  that  all  four  blade  models  give  good  correlations 
and  that  the  spring  model  slightly  overpredicts  damping.  This  overprediction,  though  not 
appreciable  quantitatively,  is  consistently  observed  throughout  the  data  range.  In  Fig.  15b 
(00  =  3°  and  =  16°)  none  of  the  blade  models  predict  the  trend  of  the  data  of  decreasing 
damping  with  increasing  advance  ratio  for  the  entire  data  range  (0  <  /:/  <  0.35).  All  the 
models  except  the  spring  model  provide  adecjuate  correlation  for  <  0.2  or  so.  Finally,  we 
come  to  Fig.  16  for  do  =  6°  and  a*  =  20°,  in  which  the  data  are  limited  to  low  advance 
ratios  (0  <  ^  <  0.15).  It  is  seen  that  the  modified  model  and  the  modal  model  provide  good 
correlations.  The  predictions  from  the  rigid  flap-lag  model  are  also  adequate.  The  trend  of 
the  data  is  predicted  by  the  spring  model,  but  it  consistently  overpredicts  damping.  This 
overprediction  is  due  to  the  overestimation  of  root-flexure  bending-torsion  coupling  by  the 
spring  model,  and  it  increases  with  increasing  blade  pitch;  for  example,  compare  Fig.  15a 
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and  15b  with  Fig.  16  at  {.l  =  0. 


5  Concluding  Remarks 

The  preceding  wake  modeling,  parametric  and  correlation  studies  show  that  dynamic 
stall  and  wake  effects  are  appreciable  on  control  inputs,  periodic  responses  and  damping 
levels  and  lead  to  the  following  findings: 

1.  The  dynamic  wake  model  with  nine  harmonics  (M  =  9)  and  with  three  radial  shape 
functions  per  harmonic  gives  essentially  converged  results  in  the  presence  of  dynamic 
stall. 

2.  The  effects  of  dynamic  stall  and  dynamic  wake  on  the  control  inputs,  periodic  responses 
and  on  the  damping  levels  increase  with  increasing  thrust  level  and  advance  ratio. 

3.  For  nearly  the  entire  database,  the  rigid  flap-lag,  modified  and  modal  models  provide 
adequate  correlation.  The  predictions  from  the  modified  and  modal  models  are  nearly 
the  same  and  run  fairly  close  to  the  predictions  from  the  rigid  flap-lag  model.  However, 
with  increasing  negative  thrust  levels,  as  is  the  case  for  combinations  of  high  advance 
ratio  and  high  shaft  angle  at  0o  =  0°,  the  modified  and  modal  models  provide  better 
correlations. 

4.  For  a  set  of  data  points  at  Oq  =  3°  and  high  shaft  angles  the  trend  of  the  data  of 
decreasing  damping  with  increasing  advance  ratio  is  not  well-predicted.  This  occurs  at 
relatively  low-thrust  conditions  and  dynamic  stall  is  not  a  major  factor.  Investigation 
is  continuing. 

5.  For  6q  =  0°  and  high  shaft  angles,  as  the  advance  ratio  increases,  the  effects  of  nega¬ 
tively  stalled  conditions  increase.  Under  such  conditions  the  predictions  show  increas¬ 
ing  sensitivity  to  changes  in  modeling  the  root-flexure-blade  assembly.  Consequently, 
the  differences  in  the  predictions  from  four  structural  models  increase  with  increasing 
negative  thrust  conditions. 


17 


6.  The  spring  model  also  provides  fairly  adequate  correlation  for  the  entire  database.  But 
it  overestimates  the  bending-torsion  couplings  of  the  root  flexure.  This  overestimation 
brings  in  increasing  quantitative  degradation  to  the  correlation  with  increasing  blade 
pitch  angle  Oq. 

7.  The  dynamic  stall  theory  is  keyed  to  the  airfoil-section  characteristics  of  lift,  drag  and 
pitching  moment  under  c[uasisteady-stall  and  dynamic-stall  conditions.  Improved  rep¬ 
resentation  of  these  characteristics  under  the  low-Reynolds-number  conditions  of  the 
experiment  merits  further  research;  this  offers  considerable  promise  in  further  improv¬ 
ing  the  correlation. 
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Figure  1;  Schematic  of  Blade  Section  Aerodynamics 


20 


0  3  6  9 

Number  of  Hannonics,  M 

Figure  3;  Oscillatory  Convergence  Trends  of  First  Lag-Mode  Damping  Levels  with  Respect 
to  Number  of  Harmonics  for  ji  =  0.3  and  Crja^  =  0.075. 


22 


(1) I  (Torsion  I)  W,  (Flap  I) 


Figure  4;  Effects  of  Aerodynamic  Modeling  on  Periodic  Responses  for  fx  =  0.3  and  Cj/ 
0.1.  (concluded  in  the  next  page) 
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Figure  7:  Control  Inputs  Vs.  Thrust  Level  from  Various  Aerodynamic 
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Figure  8:  Effects  of  Aerodynamic  Modeling  on  Lag-Damping  Predictions  from  Hover  to 
Forward  Flight  for  Crjo's  =  0.075. 
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Figure  9;  Effects  of  Aerodynamic  Modeling  on  Lag-Damping  Predictions  from  Hover  to 
Forward  Flight  for  CtIcTs  =  0.1. 
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Figure  11:  Effects  of  Aerodynamic  Modeling  on  Lag-Damping  Correlation  from  Modified 
Model  for  9o  =  0°  (concluded  in  the  next  page). 
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Figure  12:  Effects  of  Aerodynamic  Modeling  on  Lag-Damping  Correlation  from  Modified 
Model  for  =  3°. 
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Figure  13:  Effects  of  Aerodynamic  Modeling  on  Lag-Damping  Correlation  from  Modified 
Model  for  9o  =  6°. 
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Figure  14:  Effects  of  Structural  Modeling  on  Lag-Damping  Correlation  from  Dynamic  Stall 
and  Wake  Theory  for  Oq  =  0°  (concluded  in  the  next  page). 
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Figure  15:  Effects  of  Structural  Modeling  on  Lag-Damping  Correlation  from  Dynamic  Stall 
and  Wake  Theory  for  Oq  =  3°. 
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Figure  16:  Effects  of  Structurcil  Modeling  on  Lag-Damping  Correlation  from  Dynamic  Stall 
and  Wake  Theory  for  do  =  6°. 
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